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\ Abstract. We establish pointwise bounds for the Green function and consequent lin- 

earized stability for multidimensional planar relaxation shocks of general relaxation sys- 
terns whose equilibrium model is scalar, under the necessary assumption of spectral sta- 
^0 , bility. Moreover, we obtain nonlinear asymptotic behavior/sharp decay rate of per- 

\l , turbed weak shocks of general simultaneously symmetrizable relaxation systems, under 

small n perturbations with first moment in the normal direction to the front. 
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1. Introduction 

In this paper, we investigate the time-asymptotic stability of multidimensional planar 
shocks of general relaxation systems whose equilibrium model is scalar. 

Consider hyperbolic relaxation systems of general form 

d 



(1.1) 



t 



Q^{U,V) 



•-3 




u, € M}, V, ,q € W, with the condition 

(1.2) Re a{qy{u,v^{u)) < 
along a smooth equilibrium manifold defined by 

(1.3) E:={{u,vMu,v)=0} 

and T determines relaxation time. The first equation and the second r equations represent 
a conservation law for u and relaxation rate equations for v, respectively. The condition 
(|1.2|) implies that a perturbed solution eventually relaxes to the equilibrium state. 

The first equation for u can be approximated by a hyperbolic conservation law and a 
parabolic conservation law to "zeroth order" and "first order", respectively. Here, the 
corresponding order is with respect to parameter r determining the relaxation time. To 
"zeroth order", the corresponding "relaxed" scalar equation is 

d 

(1.4) ut+Y,muh,=o, 
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where fi{u) := /■'(u, and to "first order", the corresponding parabolic conservation 

laws is 

d 

(1-5) ut + Y,mn)., = {HH^^Ma:„ 

where 



-fiQv^igi - giQv^qu - {fi - fiqv^qu)qy^qu) , ifj = k 
(1.6) bi'' = { -^(^fUvHai - gUv^qu + {fu - fvqv^qu)qv^qu) 

+fvqv^{9u - aUv^qu + {fi - flqv^qu)qv^qu)) , if j 7^ A; 

is determined by the expansion of Fourier symbol(Chapmann-Enskog expansion) as in Ap- 
pendix El 

These approximation equations suggest us to investigate the existence of shock wave 
solutions. A planar relaxation shock wave is a traveling wave solution of satisfying 

(m, f )(x, t) = (n, v){xi — st), 



(1.7) 



lim {u,v){z) = {u±,v±), 

2— >±oo 



where the end states {u±,v±) satisfy v*{u±) = v± and u± is a shock solution of (|1.4|) . 

Such traveling wave solutions are known to exist for small amplitude profiles, see for 
example, |Liu[ lYoZj IMZlj . However, profiles of large amplitude may develop "subshocks" 
or jump discontinuities. We restrict here to the smooth and small-amplitude case. 

Stability of such multidimensional planar shock wave solutions has been studied for spe- 
cific models. Nonlinear stability of planar shock fronts for the 3x3 Jin-Xin model in two 
spatial dimension has been proved in [Li]. For a two-dimensional shallow river model. Ha 
and Yu proved nonlinear stability of small amplitude shocks in |HYj . Both of these analy- 
ses follow the approach introduced by Goodman in |Goj to treat the related scalar viscous 
case, based on energy estimates and conservation of mass, yielding sup-norm convergence 
to the unperturbed front with no rate. This method has since been greatly sharpened in the 
scalar viscous case, using shock-tracking and spectral (inverse Laplace-transform) methods 
and pointwise estimates on the resolvent to obtain asymptotic behavior and sharp rates 
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of decay for general scalar models; see e.g., However, up to now, no 

comparable result has been carried out for the relaxation case. 

In the present paper, generalizing the results of [HoZH IHoZ2j in the viscous case, we 
prove stability, with asymptotic behavior and sharp rates of decay, of small-amplitude 
multidimensional planar relaxation shocks of N xN general systems (jl.ip whose equilibrium 
model is scalar, under the following assumptions. 

Assumptions 1.1. 

(HO) P,c/^g G C'"+i,m > [d/2] +2. 

(HI) (i)cr{^j £,j{df^ ,dg^y{u,v)) real, semi-silultiplicity, for all ^ G W^, and (ii) 
(^a{df^,dg^Y{u,v)^ different from s. 

(H2) ^jdfi{u±)) real, distinct and different from s. 

(H3) 5Rc7(i ZU i^jidP,dg^Yiu±,v±) - (0, dqYiu±,v±)') < V(l + I^H for all ^ G 
, 6* > 0. 

(H4) The set of solutions of (11. ip forms a smooth manifold {us,V5),6 (zU ^M}. 

Let D{X,^) as in Definition 12.61 Section \2.2\ denote the Evans function associated with 
Fourier transform of the linearized operator about the wave, a function that is analytic 
in A for 5RA > —6, 6 > 0, with zeros corresponding with eigenvalue of L^. (For history and 
further discussion of the Evans function, see |AGJ1 IGZ^ IPZj and references therein.) 

Assumptions 1.2. (Strong spectral stability conditions) 

(Vl) D(-,0 no zeroes in > 0} except at = A = 0. 
{V2) {d/dX)D{0,0) / 0. 

{V3) A zero A*(0 of D{-,i) satisfies A*(0) = and 5RA*(|) < for ||| 

sufficiently small. (Existence and local uniqueness of A=k(^) are guaranteed by the Implicit 
Function Theorem and {'D2).) 

^ Though the analysis of |Go) also involves an approximate front location, this is determined by a zero 
residual mass condition convenient for energy estimates rather than considerations of asymptotic behavior, 
and involves errors of the same magnitude as the perturbation itself; see the discussions in [GMIIHoZlllHoZ2] . 
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Sometimes it is more convenient to write p.ip in the abbreviated form 

d 

(1.8) Ut + Y, A'{U)U^, = r-'QiU) 

where U = {u,vY, A^{U) = {df^,dgi)\u,v) and Q{U) = (0,g(n,f))*. 
Assumptions 1.3. 

(^1) (jl.Sp is symmetrizable in the sense that there exists symmetric, positive defi- 
nite such that A^A^ are symmetric for all j = 1,2, ...,d and A^dQ is symmetric, negative 
semidefinite. 

{A2) (Kawashima condition) There exists the operator K{dx) such that 

(1.9) mT/(e) = iKiOfio 

where K{C) is a skew- symmetric operator which is smooth and homogeneous degree one in 
^ satisfying 

d 

(1.10) ^a{\i\^A^dQ -YijK{i)A^)^ < -e\Cf for all ^ in M.'^ 

Remark 1.4. If (jl.Sp satisfies (^1) and the Genuine-Coupling condition that no eigen- 
vector of cr (^iJ2'^^^ ^j{df^ , dg^Y {u±, v±)^ lies in the kernel of dQ{u±,v±), then (^2) holds 
pCl lSKllMZSllZil IZS]. Moreover, conditions (A1)-(A2) imply (H3) pOISKlIZe]. 

Conditions (H0)-(H4), (Al)-(A2) are the standard set of hypotheses proposed by WA. 
Yong for relaxation systems [Yo], as adapted to the shock case by Mascia and Zumbrun 
[MZlj . As described in |MZ1[[MZ5] . (A1)-(A2) are satisfied for a wide variety of relaxation 
systems, in which case all of (H0)-(H4) are satisfied for sufficiently small-amplitude profiles 
under the single condition (Hl)(ii). 

Before we state our main theorem, we briefly go over the idea used by Goodman and 
Miller |GM| to give a formal qualitative description of the behavior of the linear perturbation 
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U{x,t) := U{x,t) — Lf{xi), where U{x,t) is a solution of and U{xi) is a shock wave 
solution. Linearizing ()l.ip about U{xi), we obtain the linearized perturbation equations 

d 

(1.11) Ut = LU:=- Y^iA^U)^. + {dQ)U 

3=1 

where 

(1.12) := dF^{U{xi)), dQ := dQ{U{xi)) 

depend only on the normal direction xi. 

We can give a heuristic approach to describe the behavior of the perturbation U = U — tJ . 
First, we approximate the operator e^f* by its formal spectrral projection 

(1.13) e^*(«)V(V;,C/o), 

onto the top eigenfunction of L^, assuming a perturbation expansion 

(1.14) K{^) = i^-i+~efi+--- 

of the corresponding eigenvalue A*. Next we apply the method of stationary phase to the 
inverse Fourier transform to obtain the approximation 



(1.15) U{x,t) = e^*C/o = I eVe^«-^c7odf 



(27r)'^-i 

- -U'{xi)6{x,t) 
where 5{x,t) satisfies the transverse convection-diffusion equation 
(1.16) 6t + a-V^6 = dw^iPVi5) 
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with initial data 

(1.17) 6o{x) = (^,f/o(x))i2(,^) = -{[u]-\0) J Uo{x)dxi. 

The following theorem shows that the formal linear approximation C/(x, t) ~ — 5(x, t)U'{xi) 
is in fact valid at the nonlinear level. 

Theorem 1.5. For fixed let U he a relaxation shock profile (|1.7|) satisfying (H0)-(H5), 
(A1)-(A2), (V1)-('D3), with amplitude \U+-U-\ sufficiently small. If\Uo-U\Li,\Uo-U\L2, 
\xiiUo — U)\L^, \Uq — U\j^[d/2]+3 < Co sufficiently small, then for arbitrary small a > 0, there 
holds 

(1.18) \U{x,t) - U{x, - 5(x,t))L.(,.) < CCo(l + t)-('^-i)/4-i/2+- 

for dimensions d > 2, with 5 as defined in (jl.l6p - ()1.17p and a, (3 as in (jl.l4p . Moreover, 
the above result holds with (7 = for dimensions d > 3. 

1.1. Discussion and open problems. Previous results [Lil IHYj on multidimensional 
scalar relaxation fronts rely on the structure of specific models and give only stability 
without rates or behavior. Theorem 11.51 by contrast applies to general equations, giving 
sharp decay rates and a detailed picture of asymptotic behavior. On the other hand, it 
relies on the assumption of spectral stability, which must be verified. 

We conjecture that conditions (2?1)-(X'3) might be verified by a singular perturbation 
argument like that of |PZj for the one-dimensional (system) case. Verification by this or 
other means is an important open problem. Another interesting open problem would be to 
remove the restrictive hypothesis (Hl)(ii) as discussed in [MZSj |MZ6] . 

The main obstruction to the application to the relaxation problem of the spectral tech- 
niques of |GM[|HoZltlHoZ2] is to treat the more singular high-frequency behavior associated 
with the hyperbolic nature of the equations. In the viscous case, the linearized operator 
about the wave is sectorial, generating an analytic semigroup, and high-frequency contri- 
butions are essentially negligible. In the relaxation case, the linearized operator generates 
a semigroup, and there is substantial high-frequency contribution. We conjecture that 
it is this difficulty that has so far prevented application of these techniques, despite their 
advantages of generality and detailed information on asymptotic behavior. 
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This difficulty was overcome in tlie one-dimensional analysis of |MZlj by direct calcu- 
lation/detailed asymptotic expansion. However, this would appear quite complicated to 
carry out in the multi-dimensional case. Here, we follow instead a simplified version of 
an approach suggested in |Z4| in the context of hyperbolic-parabolic systems, based on 
high-frequency energy estimates. This is quite general, and should find application to other 
problems with delicate high-frequency behavior; in particular, it applies with little modifi- 
cation to the case of relaxation equations whose equibrium models are systems, preparing 
the way for a treatment of multidimensional shock fronts in this case, another important 
open problem. We regard this method of treating high frequencies as perhaps the main new 
contribution of this paper. 

A point that should be mentioned is that the analysis of [HoZlt IHoZ2| was for arbitrary 
amplitude shock waves, whereas the present analysis is limited to the small-amplitude case. 
The reasons for this are two. First, there is a limitation already at the level of the existence 
problem, since subshocks may form for too-large amplitudes in general [Liuj . However, 
supposing existence of a sufficiently smooth shock profile, we face a technical problem in 
carrying out the energy estimates of Section H] in the presence of characteristics moving both 
to the left and to the right. This was overcome with great difficulty in the one-dimensional 
case in [" MZ5j : we do not know how or whether this is possible in multi-dimensions. 

As a final open problem, we mention the treatment for relaxation systems of existence 
and stability of relaxation shock layers in the small-relaxation time limit, analogous to the 
small-viscosity analysis of [GMWZ] in the hyperbolic-parabolic case. 

2. Preliminaries 
We begin with a series of preparatory steps, loosely following |MZ1|, [Z4]. 

2.1. Spectral resolution formulae. We derive the spectral resolution formula by proving 
that the linearized operator about the wave generates semigroup. 

Linearizing about the wave C/(xi), we have the linearized equations 

(2.1) Ut = LU:=- Y^i^'U),, + dQU, U(0) = Uq, 

j 
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where dQ = dQ[U (xi)), A> = dF^ [U {xi)). Taking the Fourier transform in the transverse 
directions x := (x2, ...,Xd)-, we reduce to a family of partial differential equations (PDE) 

d 

(2.2) Ut = := -{A^uy -i'^ijAW + dQU , U{G) = Uq 

i=2 

in {xi,t) indexed by frequency ^ G M*^"^, where U = U{xi,(^,t) denotes the Fourier trans- 
form of [/ = U{x,t) in X and " ' " denotes d/dxi. Taking the Laplace transform in t, we 
reduce to the resolvent equation 

(2.3) (A - L^)tl = Uo 

where U{xi,^, A) denotes the Laplace-Fourier transform of U = U{x,t). 

Definition 2.1. (a) The Green function G{x, t; y) associated with the linearized equations 
(|2lD is defined by 

(i) {dt — L^)G = in the distributional sense, for all t > and, 

(ii) G{x,t;y)^5{x — y) as t ^ 0. 

(b) The resolvent kernel G^^{xi,yi) associated with the resolvent equation (|2.3p is 
defined as a distributional solution of 

(2.4) iX-L^)G^^^{xi,yi) = 6{x-y). 
Formally, one can write 

(2.5) G(x,t;y) :=e^*5(x-y) 
and 

(2.6) G^^^{xi,yi) := (A - L^)-^6{xi - yi) 

Proposition 2.2 ( |MZll [Z4] ). Under assumptions (H0)~(H4), (A1)~(A2), L generates a 
semigroup |e^*| < Ce''°* on with domain T>{L) := {U : U,LU G L^}, satisfying the 
generalized spectral resolution formula, for some rj > rjo, 

-I j-q+ico p 

(2.7) G(x,t;y) = — — ^P.V. / / e'^'^+^'G^ .{xuyMdX, 

(27^^)'* L_ioo JiRd-i ''^ 



10 BONGSUK KWON AND KEVIN ZUMBRUN 

Proof. Performing an elementary energy estimate, under the assumption of symmetrizability 
of (11.811. we establish that 



(2.8) \A'^U\l2 < |A - A*rV°(-^ - A)f/|L2 

for a symmetrizer and all U € T^{L) and real A greater than some value A^,. If, in addition, 
and Q are asymptotically constant as xi — > ±00, then it is shown that L generates a 
semigroup e^* on L^, satisfying \e^^\]^2 < Ce'^* for some real ui [MZlj . This is done by 

()2.8p and a standard result of Henry. Therefore, by [Pa], P-1-, the inverse Laplace- Fourier 

Transform formula holds for L, e^*: 



(2.9) e^*/ = 7:^P.V. / / e^«-^+^*(L--A)-id|dA, 



' r}—ioo 

As a consequence, (|2.7|) holds in the distribution sense. □ 



2.2. The Evans function. Consider the homogeneous eigenvalue equation 

d 

(2.10) {X-L^)W = (X + i^^jA^ - dQ)W + {A^Wy = 

i=2 

and its constant-coefficient limits as xi ±00, 

(2.11) (L^ I - X)W = {dQ± - A^^^ - \)W - {A\W)' = 
or, alternatively, 

d 

(2.12) W = {A\)-^[dQ± -Y^i^.A^^^ _ x^w. 

i=2 

Definition 2.3. The domain of consistent splitting A is defined as the connected component 
of (A, ^) G C X M.'^^^ containing ^ = and A going to real +00, for which the coefficients 

d 

(2.13) {Ai)-'[dQ± - 5^^e,-4,,± - a) 

in (j2.12p have k eigenvalues of negative real part and N — k eigenvalues of positive real part, 
with no pure imaginary eigenvalues. 
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Lemma 2.4. Under assumptions (HO)-(Hl), (H3), 
(2.14) A c {(A,0 : 5ftA > -9\i\y{l + 

In particular, for |(A, ^)| > r > 0, arbitrary, A C {A : 5RA > —rj} x M'^^^, where r]{r) := 
Or"^ > 0. 

Proof. Noting that eigenvalues /x( A, ^) of coefficient (|2.13p relate to solutions of the dispersion- 
relation 



A(e) G a(^dQ±-^i^jAj^± 



A 

by the relation ^ = i^i, we find by Assumption (H3) that the coefficient has no pure 
imaginary eig envalues when SRA > -6'|C|V(1 + ICP) for ^ = (Ci, |), all E M, or equivalently 
3f?A > — + A straightforward homotopy argument taking A to real plus infinity 
then gives the result; see pZTl [Z3l [Z4]. □ 

Proposition 2.5. Under assumptions (HO)-(Hl), (H3), for (A,^) in the domain of con- 
sistent splitting A, there are solutions of (j2.10p 

(2.15) {(^+(xi;A,e),---,¥'^(a:i;A,|)} 
and 

(2.16) {ip~^^{xi;X,^),- ■ ■,ipJ^{xi;X,^)}, 

N = r + 1, which are locally analytic (in (A, ^) ) bases for the stable and unstable manifolds as 
xi —5- +00 and xi —oo, respectively, that is, the (unique) manifolds of solutions decaying 
exponentially as xi — > ±00. There are also solutions of (j2.10p 

(2.17) {?/;f(xi;A,f),.-.,Vfe(xi;A,0} 
and 

(2.18) {V'^+i (2^1 ; A, I) , . . . , ^+ (xi ; A, I) } 

which are locally analytic (in (A, ^)J bases for stable and unstable manifolds as xi —00 
and x\ +00, respectively, that is, manifolds of solutions blowing up exponentially as 
xi ±00 (not unique). 
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Proof. This standard result holds for general variable-coefficient systems whose coefficients 
converge exponentially as xi — > ±00 (a consequence of the gap and conjugation lemmas of 
Appendix E see [MZTl [MZ3l [Z3l [Z4] . □ 

Definition 2.6. (Evans function) For (A, ^) in the domain of consistent splitting, we define 
the Evans function as 

(2.19) D{\,0-= det {iff ,...,ip^,ip'^^-^,...,ip^)\^,=o. 

Evidently, the Evans function is locally analytic in (A, ^) in the domain of consistent 
splitting, with zeros of D{-,^) corresponding to eigenvalues of L|. We shall show in Section 
12.41 that hence D as well, extend analytically to a domain 

(A,0 G {A : 3?A > -r?} X M^-i, 77 > 0. 

2.3. Construction of the resolvent kernel. We next derive explicit representation for- 
mulae for the resolvent kernel ^. We seek a solution of form 



(2.20) 



$+(xi;A,|)A^+(yi;A,0 , > yi 
<^-{xr,X,i)N-{yi;X,i) , xi < yi 



where 

(2.21) $+(xi; A, = K(xi; i, A), • • •, ^tixr,i A)) G M^^^^ 
and 

(2.22) c^-(xi;A,|) = ((/.^^,(xi;f,A),---,(^^(xi;f,A)) eM^x(^"'=). 
Imposing the jump condition 

(2.23) (<I>+(2/i;A,0 ^-{yr,X,0)\ ^[^l \ = -{AYHvi), 



-N-{yi;X,0 

and inverting (j2.23p . we express for the resolvent kernel G^^: 
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(2.24) G^^^ixuyi)=< 



-(cI>+(xi;A,0 0)(cl>+ ci>-) \y,;X,0{A')-Hyi), 

xi > yi 

(0 <^>-ixr,X,0){^^ ^-)-\yi;X,C){A^)-Hyi), 

xi < yi 

Now, consider the dual equation of (|2.10|) . 



(2.25) 

where 



(L| - X*)W = 



L}W :={A^)*W' + {dQ* - A})W 



Lemma 2.7. For any W, W such that {L^ - \)W = and (L| - \*)W = 0, there holds 
(2.26) {W,A^W) = constant 



Proof. 



{W, A^W)' = {{A^)*W\ W) + {W, {A^W)') 

= ((A*/ - dQ* + A*^W, W) + {W, i-M + dQ- A^)W) 
= 



□ 



From (j2.26p . it fohows that if there are k independent solutions (p~[ , . . . , ip'j^ of (L^ — 
XI)W = decaying at +oo and N — k independent solutions 93^^-^, 99^ of the same 
equation decaying at —00, then there exist N — k independent solutions V'tv 
(L| — X*I)W = decaying at +00 and k independent solutions , of the same 

equation decaying at —00. More precisely, setting 



(2.27) 



*+(xi;A,e) = (V't-i(^i;^'0 ••• ^]^(xi;A,0) GM^x^^-'^) 



(2.28) 



^! {xi-X,i) = {'4>^{xi-X,i) ■■■ V'fc (3;i;A,e)) G 



t,Nxk 
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and 



(2.29) 



*(xi;A,|) = (*-(xi;A,C) *+(a;i;A,0) G 



r,NxN 



where tjj^ are the exponentially growing solutions at ±00 respectively, of (L^ — XI)W = 
as described above. We may define dual exponentially decaying and growing solutions i^}^ 
and via 



(2.30) 



(* $)^A^(* ^)^ = I 



We seek the Green function ^ in the form 



(2.31) G^|(xi,yi)= < 



$-(xi;A,OM-(A,0^+*(yi;A,|) , < yi 



where 



(2.32) M(A,0:= 



-M+(A,e) 
M-(A,|)^ 



and 



(2.33) := (^- §+). 

Note that the independence of the righthand side with respect to z is a consequence of the 
previous lemma. Thus, 



(2.34) G^_|(xi,yi) 



$+(xi;A,0 0) M(A,0 (*-(yi;A,0 o)" 
$-(xi;A,|)) M(A,0 (0 ^+{y,;X,Oy 



xi > yi 
xi < yi 
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Proposition 2.8. For (A, ^) G A, there hold 

' M+ (A, l)^]{xv. A, mi{yv. A, If for yi < < xi, 



(2.35) G^^^ixuyi)= < 



A: 

^ (A, O^J (xi ; A, O^fc (2/1 ; A, 1)* 



+ ^'Pk (3;i;A,|)(^fc (yi; A,0* /or xi < yi < o, 



M+ = (-/,0)($+ 1- 



(2.36) 
and 
(2.37) 

^ symmetric representation holds for yi > 0. 



d+ = (0,/) ($+ 



Proof. This follows exactly as in the one-dimensional case |MZlj . by (|2.32p together with 
Kramer's rule. □ 

Remark 2.9. Representation (j2.34p together with uniform exponential decay of <I>^, 
Proposition 12.51 and the fact that d± are bounded when the Evans function D := 
det($'^, <I>^) does not vanish yields uniform bounds 

\G,^^{x,y)\<Ce-'\^'y\, 

9 > 0, on the resolvent set p{L^), in particular (by assumption (PI)) for 9?A > —r], ij > 
on intermediate frequencies 1/R < |(A,^)| < -R, R> arbitrary. However, we shall not use 
this in our analysis, carrying out instead energy-based resolvent estimates for intermediate 
and high frequencies. We shall use (j2.34p only in the low frequency regime |(A,^)| << 1. 
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2.4. Low frequency bounds. We now examine in further detail behavior for small fre- 
quencies, carrying out at the same time the analytic extension of normal modes and Evans 
function beyond the region of consistent splitting. As in [HoZ2j , for our later arguments it 
will be important to enlarge the domain of ^ to complex values ^ G C*^^, and so we will do 
this at the same time. 

Lemma 2.10. Under the assumptions of Theorem for \{(^,X)\ sufficiently small, 
now taken in C'^~^ , the eigenvalue equation {L_^^ — X)W = associated with the limiting, 
constant- coefficient operator ^ has a basis of (1 + r) solutions, for m = 1, r, 

(2.38) W± = e^-(^'«>iy±(A,e), 
IJ,^,V^, analytic in A and ^, consisting of r "fast" modes 

(2.39) l^t. = lt + O{X,0, = {A'^r^St + 0(A,e), 

— (*-*' ^mY where 7^, are eigenvalues and associated right eigenvectors of dq±{A\.)^^{{), I^)* 
(equivalently, 7,^,5"^ are nonzero eigenvalues and associated right eigenvectors of {Q± — 
^Ej^i0^i)(^±)"V> ".ndl "slow'' mode 

(2.40) /.±+i(A, = + z| • a±) + -^(A + z| • ~a^f 

«i («i Y 

-\~eBu+o{\\+\i\t) 

and 

(2.41) VXi{\i):=R\^ + O{\0, 

where af = o|(±cx)), 0=*= = {a^, ■ ■ ■ ,a^)(±oo), = l//(±cx)), and h*'^ = h*j^{±oo), with 
a*, Vi , andbji^ as defined in (|A.22p . (jA.25p . and ()A.24p . Appendix\M Likewise, the adjoint 
eigenvalue equation (L^ ^ — X)*Z = has a basis of solutions 

(2.42) ll^± = e-^-(^'«>iF±(A,0, 
where 

(2.43) T/±(A,0 = T± + O(A,0 
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and 

(2.44) V;±i(A,0 = Lf* + O(A,0, 
where V is analytic in (A, are the left eigenvectors of 

iQ±-iY.^jAi){Air' 

associated with the nonzero eigenvalues —jJi^, and L*^ = (1,0^)"^. 
Proof. By the inversion of the expressions 

(2.45) A(0 = -ii ■ a* - ^'B*^ + ■■■ 

j^l Ml 

carried out in Appendix [A] for the dispersion curves near = 0, together with the funda- 
mental relation /x = i^i, we have R^^ = Vi{±oo) and 

,*,± 



nf.-.(X,h := -r U A2 + 0(A3) 



1 A*'-*- 

-± I cto*.±e\ I "11 (\ I „-c ;;± I c^o*.±c^2 
(«f 



- — (A + ^^ a± + CBTi) + -^(A + ie • a± + CBlfi) 



+0{\\ + \m 

^ , o ^11 Bn* c\ \ , ^11 \2 ^ ?t\ 



«i {a^r ' («i)^ «i 

-r^(^ ^^)' + ^(^~ • ( - 4 + 2ffe) + o(IA + III P) 

(«i)^ V {a^Y' 

(-4)(A + • a^) + 7%(A + i| • a±)^ 
fli («i J'^ 

-4l*^nl+0(|A + lllP)- 
«i 
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The corresponding adjoint computation yields = (1,0^)"^. (See, for example, the one- 
dimensional computation of [MZlj . which is sufficient to determine L"^^.) 

□ 

Normal modes. Consider again the variable-coefficient eigenvalue equations 

(2.46) (L^- - X)W = {dQ - - \)W - {A^W)' = 
and the limiting constant-coefficient equations 

(2.47) (L^_^- - X)W = {dQ± - A^- ^ - \)W - {A^^W)' = 0. 

We now relate the normal modes of p.46p to those of (12.470 . 

Lemma 2.11. (normal modes) Under the assumptions of Theorem for (A,^) € 
B{0,r), r sufficiently small now complex), there exist solutions W^{xi; X,^) of ()2.46p . 
in xi and analytic in X and ^, satisfying 

(2.48) W^{xr,X,i) = e>^^^W^{X,0 



for any k>0 and < 9 < 9, where /i^(A,^), V^{X,^) are as in the previous lemma. 

Proof. This is a direct consequence of the previous lemma and the gap lemma, Lemma 
El □ 

Proposition 2.12. There is a neighborhood of (0,0) in {X,(,) space (£, now complex) in 
which, for yi > 0, 



(2.50) 



G\,l{xi,yi) = G\^^{xi,yi) + G\^^{xi,yi), 
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where 



(2.51) Gl^^ixi,yi) = < 



, 2/1 < < XI 

- J2j V-'7 Cli'J (yi ; A, 0* , yi < xi < o 



Efc,,A^i(^'e>+(a;i;A,OV^^-(yi;A,0* , yi < < xi 

(2.52) G5^^-(xi,yi) = I Efc,, ^1(^,1)^^7(^1; A,|)Vifc (yi; A, 0* , yi<a;i<0 

[ Efejt^7fc(A,l)^7(a^i;A,|)^^(yi;A,|)* , xi < yi < 

where 



(2.53) |M+|,|d±|<Ci|D-i 

and £)(A,^) = det $ is the Evans function. 
Moreover, 



(2.54) 



, yi < < xi 



+ 



^C2(xi) Oy 
3M;r+i(a;i-2/i) 

xo(iAi+iiiy 

+0(6-^1^1-2/11) 



0(e-^|xi-2/i|) 



, yi < a;i < 

, a;i < yi < 
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(2.55) 



, yi<0 <xi 



xo{\x\ + \i\) 



0{e 



, yi <xi <0 

, xi<yi <0 



(2.56) Gl^ixi,yi) = CZ)-i(A,Oe-^^"+^^^i/'(xi)(l,0)+ 

0{\D-\X,^)\{\X\ + |||)e-^e A'r+iCA.OC^i-s/i)) 



(2.57) (^—^Gl^ixi,yi) = O (\D-\\,mX\ + |l|)e-^e M-+i(Ai){.i-j/i)^ 

Proof. By proposition 2.6, we decompose G into two parts and where G^ has no poles 
and G^ has all terms with poles, respectively. Thus, we have (j2.5ip and (j2.52p . Without 
loss of generality, we can set 

(2.58) If I = i7'(xi) + A(xi, A,Oe"^l^'il for c> 

where A(xi,A,^) = 0{\\\ + |^|) is differentiable and exponentially decaying in xi. This is 
possible since it solves (j2.10p . and is only bounded solution at zero frequency. On the other 
hand, we can check by inspection that 

(2.59) i^r+i = [(1,0) +e(yi,A,e)]e-^'-+^l^^l, 
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where Q{yi, A, ^) = 0([A| + 1^|) is differentiable and exponentially decaying in yi. By lemma 
3.2, all fast modes can be written as 

(2.60) i;- = e^^^W-ixuX,i) 

= e'^^"^ Vr (xi , 0, 0) + e^""^ ^'(xi , A, ^Vf (xi , 0, 0), 

where ^'(xi, A,^) = 0(|A| + |.^|) is differentiable and exponentially decaying in xi. Since all 
modes are fast except for one slow mode ipr+i, For yi < xi < 0, we have 

(2.61) Yl ^7^7* = ^7+ii";i + E V'zvi,"* 

= e'^^'+i^i {vr (xi , 0, 0) + ^{xi , A, i)V- (xi , 0, 0)) 
X e-i^r+iyi [(1, 0) + e(yi, A, i)] + ^ i^Ji'J* 

^C2(2;i) 0^ 
+ eM.-+i(xi-?^i) (^^(3;^^ A, 0^j~(a;i, 0, 0)(1, 0) 

+ W(xi,0,0)e(yi,A,O 

Vc2(a:i) 0^ 

+ (^e'^r+i(^i-J/i) X 0(|A| + III)) + O(e-^l^i-J'il) 
using (|2.59p and ()2.60p . For xi < yi < 0, since all (p^ are fast modes, we have, 

(2.62) ^^-(xi; A,0(^T(yi;A,|)* = 0{e-'\^^-y^\). 

j 

Thus, we have 

It is easy to check (j2.55p by differentiating (j2.54p with respect to yi. 
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For we have, for yi < < xi, 

(i,fc)^(l,r+l) 

= CD-\X,i)e-''r+iy^U'{xi){l,0) 

+0(|A| + |^~|)e-^e^^^^+i(^'«)(^'i-J'i)L'-HA,|) 
+0(e-f(^i-j/i)). 

Similarly, for 2;i < yi < and for yi < xi < 0, respectively, we have the same bounds. It 
is also easy to check (j2.57p by differentiating ()2.56p with respect to yi. 

□ 



Corollary 2.13. There is a neighborhood of (0,0) in (A,^) space (S, now complex) in which, 
for yi > 0, 



(2.63) 



GlAxi,yi) 





g Re l^~+i{xi-yi) 

xo{\x\ + \i\) 

+0(e-e|xi-s,i|^) 
C)(e-«ki-yi|) 



, yi < < xi 



, yi <xi <0 
, xi<yi <0 



(2.64) Gl^{x^,y,){0,lY = 0{\D-\X,mM + \i\)e-^ 

Proof If we multiply ([23il) by (0,/)*, then we have the result. 



□ 
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Remark 2.14. Note that G-' -(xi, yi)(0, /)* has almost the same bounds as -^G-' -(xi, yi) 
for j = 1,2. This will be important in our later derivation of pointwise Green function 
bounds, yielding that differentiation by y is roughly comparable to right-multiplication by 

2.5. Decomposition of the Green function. For fixed small (5i, ^i, r > to be chosen 
later, define a "low-frequency" part and a "high-frequency" part G^^ of G, respectively 
by 

G\x,t-y) : = 



where Fq = [—r] — ir,r] — ir] U [?] — ir, r/ + ir] U [2r] + ir, —ij — ir], and 
G'\x,t;y) := 

— -^P.V. / /_ e^^<'-y^^''G,^ixuyi)didX. 

[2mf J-ei-ioo J\i\>Sior |3A|>r 

Then, by the spectral resolution formula (|2.7p together with Cauchy's Theorem, we have 
a decomposition formula for G{x, t; y) of 

(2.65) G{x, t; y) = G\x, t; y) + G'\x, t; y). 

3. Low-frequency estimates 
We estimate the low frequency part of the Green function following |HoZ2] . 
3.1. Pointwise bounds. The stationary phase approximation 

U ^U'{xi)5{x,t) 

of (jl.l5p can be expressed alternatively (recalling that the low-frequency part of G generally 
determines large-time behavior) as 

(3.1) G\x,t-y)^U'{xi){[u\-\^)g{x-y,t), 
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where g{x,t) denotes the Green function for the constant-coefficient (d — l)-dimensional 
equation (I1.16P approximately governing normal deformation of the front. Using the analysis 
of |HoZl| together with the low- frequency description of the resolvent kernel in Proposition 
13. H we now establish the following pointwise description of the low-frequency part of the 
Green function , sharpening the formal prediction of (j3.ip . 

Proposition 3.1. Under the assumptions of Theorem for yi > 0, \a\ < 1, and some 
T], M > 0, 



(3.2) D'^G^ = X{|,,_,,|<|,+i|}X{o>,,+aM^."f^'(^i)(M-\o)5+(5,t;y) 

+ X{x,>o}^yK+ix,t;y) + R+, 



where 



d-l (S-H-<i+(Hi,t)t)*/3_^l(a-B-a+(i;i,t)t) 

(3.3) t{x, t; y) := Cgt-—e « 



and 



(3.4) K+{x, t; y) := CQit-^e s 



are {d — 1)- and d-dimensional convected heat kernels, respectively, with 



(3.5) «+:=(i_^)a + ^a^ 

1^1 ^1 1*^1 ^1 



(3.6) ^+ := (1 



\Vi\ 
\a^t 



■)/? + 



\yi\ 

\a'lt 



^11,+ 



(a^ 



a+6**)(a^ 



4b* J) ^ 



(3.7) 



B* := B* - b*b*\ 
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where a+ is given by a+ =: {af,a'^) = (d/i (n+), d/|(M+), ...d/^ («+)), b* , B* , B* , a, and (3 
are as in ()A.26P and (jC.Sp , and the (faster decaying) residual terms and 0^ satisfy 

(3.8) \Rt\<X{^.,.y,^<^att\}t-'^-'^ia + t)-'^+^ie-^^^^^^ 



[g; — y-a+tp d-1 |q 

+ Ce Mt t 2 2 e 



+ Ce -t r 2 -- ((1 + ty2X{^^>o} + e-^') 



e 



-vii\xi-yi\+t) 



A symmetric description holds for yi < 0. 

Before proving Proposition [3TT1 we first establish the following lemma, a simplified version 
of Proposition 2.8 in |HoZ2] . allowing us to vary ^ in the complex plane. Following [HoZ2] . 
denote ^ E C by 

(3.9) e = ei + ^6. 

Lemma 3.2. Under the assumptions of Theorem M.li for some 62, r?i > sufficiently small, 
\x — y\ < M{t + \xi — yi\), and any \S,2\ < S2, 

G\x,t-y) = -1^ / / e*«-(^-*)+^*G;,,(xi,yi)dAd| 



(3.10) 



Proof. By Cauchy's Theorem, it is equivalent to show 

(3.11) / / e'^-^^-y^+^'G^-Axi,yi)dXdi=0{e-'^^^\''-y\+'^). 

By assumption (2^1) and continuity, there are no zeros of the Evans function for |^i| = 5i 
and 5RA > —9, some ^ > 0. Thus, by Cauchy's Theorem again, (j3.1ip is equivalent in turn 
to 

r p—ri+ir 

(3.12) /_ _ / e^«-(^-^)+^*G;,^^-(xi,yi)dA(i| = 0(e-''i('^-^l+*)). 

J \ii\=fil,\i2\<S2 J —rj—ir 

Taking << (^1, we have on the domain of integration by our previous bounds (see for 
example Remark 12. 9p that 

\G^-^{xi,yi)\<Ce-'^\^^-y^\ 
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for 02 > 0. Taking 62 « O2, r], we thus have 

yielding the result for r/i := (1/2) inm{62, rj}. □ 
We have also the following weakened version of Proposition 2.7, |HoZ2j . 

Lemma 3.3. Under the assumptions of Theorem M.li for M sufficiently large and \x — y\ > 
M{t + \xi-yi\), 

(3-13) \G'ix,t;y)\ < C -. 

^j=2i^ + \xj-yj\) 

Proof. We first consider the simplest case of dimension d = 2. Moving ^2 from to c := 
h{x - y)/\x - y\ in 

I I e'i-^''~y^+^'G^,{x,,yi)dXdl 

62 > fixed, by the argument of Lemma [312] yields a change of order 

Jo ~ 1 + |5 - y| ' 

On the other hand. 



f I e''^-^^-y'^+^'G^Axi,yi)d\dl<Ce-^^\^-y\e^ 

J\iT\<ST.e'2=cJrn 



-S2\x-y\^Ciit+\xi-yi\) 
< ^^-{S2/2){t+\x-y\) 



since by assumption \x — y\ >> {t + \xi — yi\). This completes the proof for dimension 
d = 2. 

For dimensions d > 2, we proceed by induction, moving one component of ^2 at a time, 
starting with the component for which \xj — yj\ is largest, without loss of generality, j = 2. 
At the first step, then, this yields a change of order 

^-Vit+\xi-yi\) / g-z\x2-y2\^^ < h_ L 

Jo ~ l + |2;2-y2| 

times the maximum of a family of {d — 2)-dimensional integrals in (^3, . . . of similar 
form to the original one, plus a new integral of negligible order 0(e~^^2/^^'-*^l^~^l^). Moving 
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the next component similarly, in each of the family of {d — 2)-dimensional integrals, yields 
a factor 

C 

1 + ^3 - ysl 

times the maximum of a family of {d — 3)-dimensional integrals of similar form, plus a new 
integral of smaller order 0(e~''l^'3~^^l) times another {d — 3)-dimensional integral of similar 
form. Continuing this process, we obtain the result. □ 

Remark 3.4. For |x — y| > Mt, M sufficiently large, we have G = by finite speed 
of propagation. Thus, (j3.13p reflects a certain inefficiency in our splitting scheme. Note 
that the righthand side is time-exponentially decaying in L^, p > 1, whereas usual error 
terms 0(e~''(*"^l^~^l)) are time-exponentially decaying in all norms; thus, for practical 
purposes it is almost but not quite optimal. For the present analysis, in which we consider 
only 2 < p < oo, it is harmless. 

Proof of Proposition \3.1[ With our preparations, this follows now by exactly the argument 
used in [HoZ2j to establish the corresponding bounds on the full solution operator for 
|x — y| > Mt, M >> 1. (Note: what is actually estimated for this regime in pipZ2j is the 
low-frequency part, with the rest shown to be negligible.) 

Specifically, we note that the description of the resolvent kernel in Proposition 12.31] agrees 
with that for the viscous case in Proposition 2.5, |HoZ2j in the sense that the principal 
terms are identical, with the rest consisting of fast-decaying (0(e~^l^^~^^l)) terms leading 
to negligible errors. The only difference in the relaxation case is that there are more (at 
most r) of the latter terms than in the viscous case (at most 1). 

Then, the rest of the proof goes word for word as in the arguments of |HoZ2| . Sections 
3, 4, and 5, based on careful stationary phase estimates of the various terms. For this 
(complicated) argument, we refer to |HoZ2j . □ 

3.2. LP bounds. From Proposition 13. H we readily obtain the following bounds on 

the solution operator. 

Lemma 3.5. For t > 1 and a is a multi-index with \a\ < 1, there holds 
(3.14) \D-^G^{x,t;y)\ < C(p)t((^-i)/2)(i-VP)-hl/2 
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for all p > 1, with C ^ oo as p — > 1. Moreover, for f G and p > I, there holds 



(3.15) 

and 



[ G\x,t;y)f{y)dy <C{p){l + t 

J LP{x) 



-({d-l)/2)(l-l/p)[^| 



(3.16) 



Gl^{x,t;y)f{y)dy < C(p)(l + t)-«'^-i)/2)(i-VP)-i/2 [_^|^^ . 
LP{x) 



Proof. By calculating the norm of the (d— l)-dimensional heat kernel, together with the 
expression of (|3.2p . we have 



(3.17) 



< (j^i{d-i)/2)ii-i/p)-\a\/2 ^ C7(p)e-^*, 



r] > 0. The inequalities (j3.15p and (j3.16p follows from (j3.14p and the triangle inequality. □ 



Lemma 3.6. For p > 1 and C{p) as in Lemma Vj. 5\ 



(3.18) / G\x,t,y){Q,Ir)'f{y)dy < C(p)(l + t)-('^-i)/2((i-i/p)-i/2) |;|^, . 

J LP 

Proof. By (12.631) and (I2.64p . we know that ^(xi, yi)(0, /)* has the same bound as ^G'^ ^(^^i) yi)- 



Thus, 



fG'{x,t,y){0,Iryfiy)dy 



has the same bound as 



LP 



fGUx,t;y)f{y)dy 



()3.16p . we have the result. 



LP{x) 



■ By 

□ 



Lemma 3.7. The low frequency part of Green function G^ {x, t; y) associated with {dt — L) 
satisfies 



(3.19) 



[G'{x, t; y) - U'{xr)Ii\xi)g{x - y, t))f{y)dy 



LP{x) 



for all t > 0, where n(xi) is the left zero eigenvector dual to the right zero eigenvector 
U'{xi) ati=0 
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Proof. By ()3.2p . we have 
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(3.20) 



\G\x, t; 0, y) - U\xi)U\xi)g{x - y, t)) f{y)dy\L,^,) 

< Ct-{(rf-l)/2){l-l/p)-l/2_ 



Then using (13.20p . we have 



(3.21) 



LP{x) 



G'{x,t;0,y) - U' {xi)n\xi)g{x - y,t))f{y)dy 
<\ j G'ix,t;0,y)-U'ixi)n\xi)gix-y,t))fiy)dy\LP^^)\f\Li(.) 

< (j^-{{d-l)/2){l-l/p)-l/2 ^ 



On the other hand, we have 



(3.22) 



LP{x) 



G'{x,t;yuy) - G'{x,t-0,y))f{y)dy 

= j j 9y^G^{x,t;9yi,y)yif{y)d9dy 

< sup \dy^G^{x,t;y)\LP(x)\xif{x)\LP{x) 
y 



LP{x) 



Combining (j3.2ip and (j3.22p together with the triangle inequahty, we have the result. □ 

4. Damping and high frequency estimates 

We now carry out the main new work of the paper, establishing high-derivative and 
high-frequency bounds by energy estimates following the approach introduced in |Z4] for 
the hyperbolic-parabolic case. 

We denote U = {u, f )* for our convenience. Let U he a solution of (jl.ip and f/ be a 
traveling wave solution of (jl.ip . Define the nonlinear perturbation U{x,t) := U{x,t) — 
U{x, t) where fj{x, t) = U{xi - 5{i, t)). We also denote A = A{tj), A = A{U), dQ = dQ{U) 
and dQ = dQ{tJ). 
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Lemma 4.1. Multi- dimensional scalar relaxation equations (jl.ip can be put in a quasilinear 
form as follows: 

(4.1) Ut + Y,^'U,^-dQU 

j 

= {dt - L)6{i, t)U'{xi) - R^, + (0, IrfNiU, U)-M=:f 

where 

d 

(4.2) M = Y,{A' - Ai)U^^ = 0{\U\\U'\\Vi5\), 

i=i 

(4.3) R^, = 0{\5t\\U'\\5\ + |Vi.5||C7'I|5|),, = 0(|5|(|5t| + \Vs^5m"\) 
and 

(4.4) N{U,U) = 0{\U\^). 

Proof. Consider the multi-dimensional scalar relaxation equations. 



d 



P{n,v)\ / 



(4-5) +E 



For [/, there holds 
(4.6) Ut + Y, 



P{U)\ i o 



=0 



-idt-L)6{x,t)U'ixi)-R,„ 
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where 



(4-7) R^.=-^tu'\:\-s+Y.(-^^.)[', 



U' 



XI 

XI —5 



The last equaUty is true since the stationary shock wave solution U'{xi) satisfies 
'f\U)\ ( 



Xl 







(4.8) = \ _ U'{x,) - 1 t7'(xi) 



dg\U{xi))j j \dq{U{xi)) 



The last equation is obtained by multiplying the shock wave equaition by 6 and adding 
to the second last equation. 
Let U he a solution of 



(4.9) Ut + Y, 



1 



PiU)\ ^ I 



After subtracting ()4.6p from (j4.9p . we have the perturbation equation for U 

d 

(4.10) Ut + Y, iA'{U)Ux, - A\^)U.,.) - {Q{U) - Q{U)) 

i=i 

d _ 

= Ut + Y, {F^iU) - F\V))^^ - {Q{U) - Q{U)) 
i=i 

= {dt - L)6U' - R^,. 
Keeping a quasilinear form on the left hand side, we have 

(4.11) Ut + Y,^'U,^-dQU 

j 

d 

= {dt - L)5{x, t)U'{xi) - R,, + (0, IrYN{U, U) - ^(i^' - 1^)% 

i=i 

= {dt - L)6ix, t)U'{xi) - R,, + (0, IrfNiU, U)-M=:f 
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where 

d 

(4.12) M = Y,{A'-^')U,^ = 0{\U\\U'\\Vi6\), 

i=i 

(4.13) i?,, = OiMU'\\5\ + \V^S\\U'\\5\).,, = OmiStl + \V,5\)\U"\) 
and 

(4.14) N{U,U) =0{\U\^) 

since A^-A^ = A^{U)-A^{U) = /J dA^ [U +e{U -U))Ude = 0{U) and U^^ = -U'6^.. □ 
Assumption 4.2. T/ie operator K{dx) is defined by 

(4.15) mT/(e) = ii^(o7(e) 

where K{(,) is a skew- symmetric operator which are smooth 
and homogeneous degree one in satisfying 

d 

(4.16) ^a{\i\^A^dQ - ^ijK{i)A^)^ < -e\C\^ for all ^ in M.'^ 

i=i 

Remark 4.3. This is a standard assumption of Kawashima which is satisfied when (jl.Sp is 
simultaneously symmetrizable and satisfies genuine-coupling condition. (See |Yo| and |Ze|.) 



Proposition 4.4 (Damping estimate). If |[/|//s(t) < e sufficiently small for < t < T 
where s > [^] + 2, then, for some 6 > 0, there holds 

(4.17) \U\ls{t) <e-''\U\l.iO) + C f e-'^'''\\U\l,{s) + €{s))ds 

Jo 

for < t < T , where 

(4.18) e{t) < C\\/t,^6\l2 = CCiil + t)-'-^-\ 



Proof. Let a be a multi- index with \a\ = r > 1. Taking a differential operator (9" on the 
equation (|4.11|) . we have 

d 

(4.19) d'^Ut + J2 {A^U,,) - (dQU) = d^f 
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where / := {dt - L)S{x, t)U'{xi) - R^^ + (0, IrYN{U, U) - M. 

Taking the L inner product of A^d^U against d^U, we have the energy estimate: 



< lii^t + ^iA'A^h,)d^U,d^U) + {iA^dQ)d^U,d^U) 

j 

k=l j 

\a\ 

+ J2 {A\d',dQ)idt\-'U), d^U) + {A'd^f, dPU) 

k=l 

< C{\U\wi,oo + \U\\^r,^) \U\jjr + e\U\jjr + C\U\l, 

+ {A^dQd^U,d^U) + eit), 

< C(|[/|h.1,co + |i7|^^l.oo) \U\%r + e\U\jjr + C\U\l2 

+ {iA''dQ).d^U,d^U) + eit), 



where 



eit) = 0{\Vt,^6\l,) = C{l + t)-^-\ 

The second last inequahty is true by Moser's inequahties and Sobolev inequahties. For each 
a with |q:| = r > 1, we define a := a — ej where j = min { i : cij is maximal}. Then, 
|a| = r — 1. Let Sr '■= {{a, a) : \a\ = r}. Similarly, taking the inner product of 
K{d^~°')d^U against d^U , we have the auxiliary energy estimate: 

(4.21) \j^{K{dr^)d^U,d^U) = ^±{^Kie-''){iO''U,{iO''U) 

= {iK{e-^){iO^USifUt) 

3 

using Plancherel's inequality together with the equation 

(4.22) u^ = -Y,iCjAiU + H, 
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where 

(4.23) H ■.= Y,i^-- ^^)U+idQ)U + f. 

j 

By a direct calculation with the Moser inequality and the assumption of smallness of U, we 
have 

(4.24) |(^0°^|L2(c) = \d^H\L2^.) < C\U\L^\U\Hr-. + e\U\Hr-i. 
Thus, we have 

\a\=r 

d 

{a,a)GSr i=l 

+ e\U\lr+C\U\l,+e{t) 
<-OYl (^'"^' ^> + + C\U\l + e{t) 

\a\=r 

< -e\U\lr + e\U\lr + C\U\l2 + e{t) 
The second last inequality is true by (|4.16|) . By (|4.25|) and choosing e < 9/2, we obtain 

(4.26) 1±(^{A'U,U) + j2 E c-'{{A'd^U,d^U) + {K{dr^)d^U,d^U))) 

r=l \a\=r 

s d 

^E E c-^(e'"(e'("-")(^°Q)--E^^-^(^""")^-)^'^> 

r=l (a,a)G5r j=l 
r=l 

<-^,\U\ls+C\U\l,+e{t) 
so long as |?7|vi/i.oo is small. We define 

s 

(4.27) 8{t) := {A'U,U) + Y, E c'"- {{A^d^U,d^U) + {K{dr^)d^U,d^U)) 

r=l \a\=r 
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It is easy to check that £{t) is equivalent to \U\jjsit). Using i^^l6\i and (g^ZI), we have the 
Gronwall-type inequaUty 

(4.28) j^£{t) < -e£it) + C{\U\Ut)+e{t)). 

Therefore, we have the result. □ 

Lemma 4.5. (High-Frequency Operator Estimate) Let G^^ be the high-frequency part of 
Green function associated with {dt — L). For any f G there holds 

(4.29) I j G'\x,t-y)f{y)dy\^^ < Ce-^Vlffs. 
for some 6 > 0. 

Proof. First, we establish the "high-frequency" resolvent bound by Kawashima-type energy 
estimate. To this end, consider the eigenvalue equation 

(4.30) (A - L^)U = f 
Set 

(4.31) (^°(A - L^){1 + 111 + dM (1 + III + d^i)U) 

= {A\l + III + d,,)f, (1 + III + d,,)U). 

We obtain by the simpler linear version of Kawashima-type energy estimate as in Proposition 
14.41 (the damping estimate) , that 

(4.32) Re A|(l + ||| + d,,)U\l,^^^^ < -0|(1 + ||| + d,,)U\l,^^^^^ 

+ C%l + \i\ + d^^,)f\l^^^^ + G*\U\l^,^^y 

Taking the imaginary part of the inner product of U against XU = L^U + /, we have 

(4.33) |ImA||C/|i.(,^)< 

C2\f\l^,,) + (2e)- Vli2(,,) + e|(l + III + 5.J^|i.(,^). 
By (j4.32p and (j4.33p . we establish the high-frequency resolvent bound 

(4.34) I (A - i^-)''|j,i(,,) < C for |(|, A)| > R and mX > -9, 
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for some R,C > sufficiently large and 6 > sufficiently small. 
Moreover, we have the intermediate frequency bound, 

(4.35) [(A - L^)-^\hi{xi) < C for R-^ < |(|, X)\ < R and mX > -9, 

for any R > and C = C{R) > sufficiently large and > sufficiently small. This follows 
by compactness of the set of frequencies under consideration together with the fact that the 
resolvent (A — is analytic with respect to in A). 

On the other hand, it is easy to check the following resolvent identity using analyticity 
on the resolvent set p{L) of the resolvent (A — L^)~^, for all / G P(L|), 

(4.36) (A - L^r'f = X'HX - L^r'Lp + X'^L^f + X'^f 

Using the resolvent identity ()4.34p . the high-frequency solution operator S^^ can be written 

as 



/—di+ioo p 



i\ — ioo 



xe 



{X-L,y'f{xi,OdidX 



P.V. 



+P.V. 



„ . L_i ^l?P+|aA|2>ei+e2 

■Ol—ioo JK" ^ 

xe^5 (*-j/)+^*A-2(A - L^)-^Lp\xi,i)d^dX 

— 01+ioo 



X 



e'^<^-y^+^'X-^Lj{xi,i)didX 



/— 6i+joo /• 

xe'^-^^-y'^+^'X-^fixuOdidX 
=: A + B + C. 

For f e H^, there holds 

(4.37) l^l^i < Ce-^i* sup ,|(A - L^-)-^|^i |/|^3 



< C*e-^''\f\ 



H3. 
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Similarly, we have 
(4.38) 



Using the triangle inequality, we have 

-— 6i+ioo 



(4.39) 



\C\m < 



P.V. / A-^e^^dA / e"-V(xi,Ode 

—di — ioo 



+ 



' —Oi—zr 

<2r/eie-''y\Hi. 



Hi 



P.V. / X-^e^^dX = I 

2Tn J_ei-joo 



The last inequality is true since the first term is zero, by the inverse Laplace transform 
identity 

6it>^ 

1 Oit < 0. 

By the expression of S^^ , together with (I4.37p . (14.380 and (|4.39p . we have 
(4.40) \S''\HS^L4t)<Ce-''. 

Therefore, we have the result. □ 

Remark 4.6. The argument of Lemma [4.5l based on the initial-value problem with homo- 
geneous forcing, greatly simplifies the treatment in |Z4] based on the initial-value problem 
with homogeneous data and inhomogeneous forcing. The argument is quite general, in par- 
ticular extending without modification to the case that the equilibrium model is a system. 

5. Nonlinear decay and asymptotic behavior. 

We now carry out the nonlinear analysis by an argument combining the approach of 
|HoZ2| in the low-frequency domain with that of |Z4] in the high-frequency domain. Let 
the approximate shock deformation 5{x,t) be a solution of the constant coefficient equation 



(5.1) 

with initial data 
(5.2) 



5t + d-Vs;6 = div£(/3Vi(5) 



6o{x) = -i[u]-\0) / U{xi,x,0) -U{xi)dxi. 
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Here, a and /? are as in Appendix lAl 

We define a smooth approximation of 5, denoted by 6^, by 

(5.3) 6%;t):=v'*S{;t) 

where rj^ is a smooth mollifier supported in -6(0, e). Note that 5^ satisfies the same convected 
heat equation as 5 does with C°° initial data 6q = t]^ * 6o. Define 

(5.4) U{x, t) := U{x, t) - U{xi - 6%x, t)). 
Lemma 5.1. For [Uq — U\li < Co; multi-index \a\ < K , there holds 

(5.5) |a?5=(-,t)|i2 < C7Co(l + t)-('^-^)/^-l"l/2 
and 

(5.6) |5^(-,t)-5(-,t)|z.2 <CCoi-("-^)/'-'/' 
where C = C{e,K) is a constant. 

Proof. Letting g{x,t) be a Green function for (j5.2p . we have 

(5.7) 5%;t)=g{;t)*6%;0) 

and 

(5.8) 6%;t)-d{;t)=gi;t)*i5%;0)-5i;0)). 

Thus, we have 

(5.9) m6^;t)\L2<\d^g{.,t)\L2\6'{;t)\Ll 

<\dig{;t)y\v'\LAS!,\L^ 
by the standard fact that a heat kernel g decays as 

(5.10) \d2g{-,t)\L2 < Ct-^d-l)/4-\a\/2 
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and observing that \6q\ii < C\IJq — tJ\ii < CCo by definition of Sq. On tlie otlier liand, 
liave also 



(5.11) 



\diS%;t)\r^2<\gi.,t)\L,\d^rj%2\6i;t)\L. 

<CCoe-l"l. 



By ()5.9p and (|5.1ip . we have the first claim. Expressing 



(5.12) 



6%;t)-6i;t) = ig'i;t)-gi.,t))*6i;0), 



and noting that 

(5.13) \ig'{y,t)-giy,t))\L2^y)= j {g{y - ~z,t) - g{y,tM{~z)d~z 

-1 

{yx9{v-0z,t)-z)T]'{z)dedz 



< / \^xg{y - 0z,t)\L2^y)\zr]''{z)\Li^-^)de 
Jo 



In the last inequality, we have used the fact that \zrj^{z)\ii^^'j < Ce. Thus, we obtain 



(5.14) 



\5'{;t)-5{;t)\L2 < \g'{;t)-gi;t)\L2\6{;0)\Ll 



Lemma 5.2. // \Uo - U\li, \Uo - U\l^, \xi{Uo - U)\l^ < Co- Then, 

(5.15) |C/o|ii, |f/o|L-, Ixif/ol^i < C{e)Co. 
Proof. We have 

(5.16) C7(xi) - U{xi - <5§) = / U'{xi - e5l)5lde. 

Jo 



< 







40 BONGSUK KWON AND KEVIN ZUMBRUN 

Then we have 

(5.17) \u(^:,^)-U{xi-6'o)\LPi.) 

Jo 

-1 

\U'\ LP {xi)\^o\ LP {x)d(^ < CCo 

By (j5.17p together with the triangle inequality, we have 

(5.18) \Uo\lp(x) < \Uo - U\LPix) + \U{xi) - U{xi - 5^)|lp{:.) < {C + l)Co. 
Similarly, we can prove the other inequality. 

Lemma 5.3. The nonlinear residual U{x,t) defined in (|5.4|) satisfies 

d 

(5.19) Ut + Yl (^'(^)^)., - 

d d 

where 

(5.20) R,, = 0{\6l\\U'W\ + \V^5'\\U'W\).,, = 0{Wm\ + \Vx6'\)\U"\) 
and 

(5.21) = 0{\U\^) and = 0{\5^\\U'\\U\) for j = G,l,...,d. 
Proof. Let f7 be a solution of 

d 

(5.22) Ut + Y. ^'(^)U^, - Q{U) = 0. 
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For U{x,t) = U{xi — 5^{x,t)), there holds 
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_ d. _ _ 
(5.23) Ut + Y,A'{U)U.,^-Q{U) 

i=i 

d 

= -5llJ' + Y,A\U){-6l^)U' 

d 

= -5lU'{x^) + Y^A^{U{xM-5l,)U'{x{) + R 
i=2 

= -{dt-L)5'tJ' + R^,. 



If we subtract (j5.23|) from (|5.22|) . we have 



d 

(5.24) Ut + Y, {AHU)Ux, - A\U)%.) - {Q{U) - Q{U)) 

d _ 

= Ut + Y, {F^iU) - F^{U))^^ - {Q{U) - Q0)) 
i=i 

= {dt - L)6U' - R^, 



By Taylor expansion of about U and (j5.24p . we have 



d _ 

(5.25) Ut + Y, {^'iU)U)^^ - dQ{U)U 

d 

= {dt - L)6U' -R^,+Y % + (0' IrfN'^iU, U) 
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By Taylor expansion of about U and (|5.25p . we have 
d 

(5.26) Ut + J2 {^'iU)U),^^ - dQ{U)U 

d 

= {dt - L)6'U' -R.,+Y, ^4 + (0' IrYN' 

i=i 

d _ 

+ ^ ({A^{U) - A^{U))U) - {dQ{U) - dQ{U))U 

d d 

= {dt - L)S'U' -R.,+Y^ % + (0, IrYN' + Si, + (0, 

i=i i=i 

where 

(5.27) R^, = om\u'\m + \v,nu'\\5'\y,, = ommi + iv^s^iwd, 

(5.28) = 0(lC/p) for j = 0, 1, d 
and 

(5.29) = 0{\6'\\U'\\U\) for j = 0,1,..., d. 



Lemma 5.4. For f £ C D L'^ , there holds 

(5.30) /* / G{x, t - s; y){ds - Ly)f{y, s)dyds 



fix,t)- / G{x,t;y)fiy,0)dy 



Proof. Integrating by parts, we have 

r-t-e 







G{x,t - s;y){ds - Ly)f{y,s)dyds 
G{x, e; y)f{y, t - e)dy - j G{x, t; y)f{y, 0)dy 

(5.31) + V^^' ~ * " '^"^y"^'- 

By duahty and letting e — > 0, we obtain the result. 



□ 
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We are now ready to prove our main result. 
Proof of Theorem \1.5[ By Lemma 15.11 it is sufficient to show 

(5.32) \U{x,t)-U{xi-6^x,t))\^,^^^ < C7Co(l + t)-("-^)/'-^/'+^ 
Define 

(5.33) C(t) := sup o<s<t\U{; 5)|l^(1 + s)('^-i)/4+i/2-<.. 
Claim 5.5. For all t > 0, there holds 

(5.34) at) < Ci{e{t) + Com + cd < c2(Co + c\t)). 

From this result, it follows by continuous induction that 

(5.35) C{t) < 2C2C0 for all t > 0, 

provided Co < i-e., C(*) remains small for all t > 0. 

Proof of claim. Applying Duhamel's principle, we can express 

(5.36) U{x, t) = {j G{x, t; y)Uo{y)dy - J G{x, t - s; y){ds - L)6'U'{y, s)dyds 

+ f j G'{x,t - s;y){-Ry, + Y,N'y, + Y.%^(y^'^^y'^' 

Jo J j=l j=l 

+ f I G''{x,t - s;y){-Ry, + Y,W^^ + ^Sl^){y,s)dyds 

+ 1^ I G'{x, t-s-y) ((0, 1,)*iV° + (0, IrfS"") (y, s)dyds 

+ 1^ I G'\x, t - s; y) ((0, /.)*A^° + (0, /r)*5°) {y, s)dyds 
= I + II + III + IV + V. 
First, we establish that 

(5.37) \I\l2 < CCo(l + t)-('^-i)/4-i/2_ 
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Using ()3.19p . we have 

G\x,t-y)U^{y)dy 



(5.38) 



L2 



< 



tj'ixi) I g{i-y,t){U' I U{yi,y,0)dyi)dy ^^(x) 



= U'{xi){6'{i,t)-6{x,t)) 



L2(x) 



+ ccot- 



-l)/4-l/2 



By standard semigroup theory, we obtain the short-time bound 



(5.39) 

By ([OS]) and (|09]) . we have 
(5.40) 
By (jl29D, 



G{x,t;y)Uo{y)dy 



< C\Uo\l2 < CCo for t < 1. 



G'{x,t;y)Uo{y)dy < CCo{l + t)-^''-'^/'-'/' . 



(5.41) I I G''{x,t;y)Uo{y)dy\^^ < e-''\Uo\HS < GCo{l + t)-^''-''^/''-'/^ 



On the other hand, by (15.30p . we have 



(5.42) 




G{x, t - s- y){ds - L)5'U'{y, s)dyds 



L2 



L2 



5'U'ix,t)- J Gix,t;yWU'){y,0)dy 

< 6'U'{x,t)-U'{xi) I g{x-y,t)6%y,0)dy 
= GCot-^''-'^^'-'/'. 

Combining this with the short time bound, we obtain 



L2 



+ CCof 



-(d-l)/4-l/2 



(5.43) 




G(x, t - s- y){ds - L)5'U'{y, s)dyds 



L2 



< CCo(l + t) 



-(d-l)/4-l/2 
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We now establish that 
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(5.44) 



\II\LP<CCo{l + t) 



-{(d-l)/2){l-l/p)-l/2 



Using (|3.16p . (|5.5p . and (|5.2ip . together with the definition of C(i), we have 
(5.45) \IIa\L^ ^\Io I ^-.'^'(^'^ - s;y)N^{y,s)dyds 



L2 



<C / (l+t-s)-('^-l)/4-l/2|C/j2,(s)ci5 

Jo 

< Ce{t) [\l + t - s)-('i-l)/4-l/2(l + ^)~(rf-l)/2-l+2a^^ 
JO 



and 
(5.46) 



L2 



|//fe|i2 = j j d:r,G\x,t - s;y)S^{y,s)dyds 

<C !\l+t-s)-^''-^^/^-^/^\S^L^ds 
Jo 

< cco Ai + 1 - sr'^''-'y'-'/'\6'\L2{s)\u\L2{s)ds 

Jo 



The last inequality is true due to the following calculations; 

rt/2 



(5.47) 



(1 + t _ ,)-(rf-l)/4-l/2(^ ^ s)-{'^-l)/4(l + 5)-('^-l)/4-l/2+<Xrf^ 

< (1 + t/2)-('^-l)/4-l/2 + ^)-(d-l)/2-l/2+.^^ 

< (l + t/2)-('='-i)/4-^/2|(l + t/2)-('^-^)/2+V2+- _i| 

< C7(l + t)-(rf-l)/4-l/2+<x|(i + ^)-{d-l)/2+l/2 _ ^ 

< C7(l + t)-('^-i)/4-i/2+<x 



It is true with arbitrary small a > for d = 2 and o" = for d > 3. 
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(5.48) [\l+t- s)~('^-l)/4-l/2(l + c,)-(<i-l)/4(l + ,)-{rf-l)/4-l/2+a^g 
Jt/2 

< (1 + i)-(d-l)/4-l/2+<x(;L ^ i)-{rf-l)/4 f (l+t- s)-(^-l)/4-l/2rfs 

< (1 + i)-(<i-l)/4-l/2+.(-L ^ i)-{rf-l)/4 f {l + t- 

Ji/2 

< (1 + ^)-(<i-l)/4-l/2+<7^-^ _^ ^^-{d-l)/2+l/2+e 

< (1 _^^)-(d-l)/4-l/2+<7^ 

The last inequality is true if we choose arbitrary small e > for d = 3 so that —{d — 
1/2 + e ^ —1 and choose e = otherwise. 

Similarly, using ([3J6]) . ([531) . and (fOOl) . 

(5.49) |//c|l2= / I d^,G\x,t-s-y)R{y,s)dyds 

< c Al + t - s)-^^-^^'^-^'^\R\l^ {s)ds 

Jo 

<c l\i+t-s)-^'-'y'-'/'\6'\MsmM') 

Jo 

< CCl [\l + t - c,)-('i-l)/4-l/2(i + ^)-(d-l)/2-l/2^^ 

Jo 

<CCo'(l + t)-('^-l)/4-V2. 

By ([535]), ([Oe]) and (fOO]) . we have 

(5.50) 1/71^2 < C{e{t) + Com + Co')(l + t)-{'^-i)/4-i/2+^ 

We now establish that 

(5.51) < CiCi + Com + CHmi + t)-('^-i)/4-V2+.. 
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By ([09]) . gUD, and dOOjl - dOTT) . we have 



L2 



(5.52) I 1^ I G'\x, t - s; y){N^^ + S^^ + Ry,)dyds 

Jo 

<C [ e-^^^-'\\U\L-o\U\Hi + \6'\l--\U\h^ + \6^\L--\Vi6^\Hi)ds 
Jo 

< cici + Coc(t) + c'mi + i)-('i-i)/4-i/2+. 

as long as and thus |C/|p^/i,oo remams sufficiently small. We shall verify in a moment 

that it indeed remains small. 

Similarly, we can establish the estimates for IV and V. By the same proof as the one for 
//, together with p.l8p . we have 

(5.53) 1/^1^2 < C{e{t) + CoC(t) + Co')(l + t)-('^~i)/4-i/2+<x. 

By an identical calculation as for the /// term except for the fact that V has one less 
derivatives, we can get 

(5.54) \V\L2it) < C(C'(t) +CoC(t) + Co')(l + 0"^'"'^/'"'/'+^ 
Therefore, 

(5.55) \U\L.it) < C7i(C2(t) +CoC(t) +Co')(l + 
We have the desired inequality, 

(5.56) at) < Ci(C'(t) + CoC(t) + Co) < C2(Co + C\t)) 
so long as \U\h^ < e sufficiently small. 

To complete the proof, we show that (i) < £ remains small for all t > if we choose 
\U\h''{0) < Co small enough, by the continuation argument. By local well-posedness, for 
sufficiently small e > 0, we can define 

(5.57) T := sup {r > : |t/|//'^(T) < e} > 0. 
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As shown above, we have 

(5.58) C{t) < C2(Co + C\t)) for t G [0, T), 
which imphes 

(5.59) C(i) < 2C2C0 for t G [0,r). 

Using ()4.17p again, we have, by choosing so smah that there holds 

(5.60) |c/|^.(T)<e-^% + ^<e. 

By continuity of i7^-norm, there exists h > such that |C/|//s(t) < e for t G [0, T + h), 
which contradicts to the definition of T. So, T = 00, i.e., |[/|//s(t) < e for t G [0, 00). Thus, 
the desired inequahty (|5.56|) is true for all t > 0, which implies that 

(5.61) Cit) < 2C2C0 for all t G [0, 00). 

This completes the proof. □ 
Theorem 5.6. If \Uq — U\ff[d/2]+5 < Co sufficiently small, then there holds 

(5.62) \U{x,t)-U{xi-6^x,t))\^^^^^ < CCo(l + t)-("-^)/'-'/^ 

Proof. Using the expression (j5.36p similarly as in the previous theorem, we get the 
bounds for each term. 

(5.63) |/|lco < C7Co(l + t)-('^-i)/2-i/2_ 
Using ([XTUj) . ([^ and PHH]) . we have 

|//a|L- < [' \Gi\Lo^{t-s)\U\l,{s)ds 

Jo 

< CCl [\l + t - s)-('^-l)/2-l/2(l + ,)-(.i-l)/2-l+2.^^ 

Jo 

<CCo'(l + t)-('='-l)/2-l/2^ 
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and 



\Ih\L^ < f\G' \L^{t-s)\U\L2{sW\L2{s)ds 

Jo 

< CCo Al + t - s)-('^-l)/2-l/2(i + ^)-(d-l)/4-l/2+<x(^ ^ 
JO 

< [\l + t - s)-('i-l)/2-l/2(l + >,)-(d-l)/4-l/2(i ^ 5)-('i-l)/4^^ 

<C7Co'(l + t)-{'^-l)/2-l/2_ 



We establish that 



(5.64) 

By (jl29]), we have 



\III\l^ <C7Co(l+t) 



-{rf-l)/2-l/2+o- 



G''{x,t-s;y)N^4yds 



< 

'0 
JO 



G''{x,t-s;y)N^4yds 



Hld/2] + l 



-e{t-s 



Xj l//[d/2]+4 



ds 



< C / e~'^^'-'^U\Lo^\U\HW2H5ds 
Jo 



<ceo{i+t) 



-(d-l)/2-l+2o- 



G^\x,t- s;y)Sl dyds 



< 



G^\x,t- s;y)Sl dyds 



jj[d/2] + l 



JO 

JO 
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j j G^\x,t - s;y)Ry^dyds 



j j G^\x,t - s;y)Ry^dyds 
Jo 

<cci{i + t)-'-''-^y^-'>/\ 



Thus we have estabhshed (j5.64p . The identical calculation as in the estimation of /// gives 
the desired bound for V, which is 



(5.65) 

Therefore, we have 
(5.66) 



I^U- <ccl{i + t) 



-(d-l)/2-l/2+cr 



Appendix A. Expansion of the Fourier symbol 
Low frequency expansion. We carry out the expansion of -P(^) in ^ about zero. 



□ 



(A.l) P{0 = dQ-iY,^j^' 







d I fj fj 
Ju Jv 



EM" 



j=l \Qu Qv J j=i \gu Qv ^ 

Claim A.l. To the second order, dispersion relations 

d 



(A.2) 

are given by 
(A.3) 
and 
(A.4) 



X{^) = a{dQ-t^(,A^), A(0) = 



X{^) = -i^.a*-eB*C + 



j,k=i 
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with a = (ai, 02, aa) and B* 
(A.5) 



j,k=l 



where 



~ fu fviv 1ui 



(A.6) 



% = < 



-\{fiQ.v^{9iu - giqv^qu + {fu - fUv^ Qu)qv^ Qu) 



(A.7) 

and 

(A.8) 

Proof. Set 
(A.9) 



2v qu J 



-qv ^qu + nv ^{ai - aiqv ^qu + a*q^ 



d d d 

j=l j,k=l 



and let =1 I for m = 0, 1, 2. Cohecting the 0th order term, we have 



(A.IO) 

which yields 



dgyo = 0. 



(A.ll) = 
Collecting the 1st order terms, we have 



-qv ^quT^ 
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(A.12) dQV^ - i{A^ - a*I)V^ = for j = 1, 2, d. 

Examining the first coordinate, we have 
(A.13) = i{fi-fiqv\u-a*iy 

So, r° is the right eigenvector of fi — fiq^^Qu corresponding to the eigenvalue a*y Let f be 
the counterpart left-eigenvector. Then, the eigenvalue a*^ is given by 

(A.14) a* = fUi - fiq-'quy. 

Examining the second coordinate equation, we have 

(A.15) Qur] + Qys] = i{gi - giq~^qu + a*q~^qu)r°, 

which yields 

(A.16) s] = -q'^quTj + iq^^igi - giqv^qu + ajqv^qu)r^- 

so. . l^'] -( 

\SjJ \-qv ^quTj + iqy ^ {gL - giqv ^qu + af-q^ ^qu)r^ j 

Collecting 2nd order terms (^j^k term), 
(A.17) + i{a*I - A^)V,' + dQV^ = 

For j = k, the first coordinate equation yields 

(A.18) b*jr^ + i((a* - fi)r] + fi{q-\ur] - iq-\gi - giq-\u + a^q-^uY'^)) 

= + ^('^j - /« + fiqv^qu)r] + fiqv^igi - dqv^qu + ajqv^quV^ = o 

For j k, we have 

(A.19) b*,V' + i{a*I - A^)V,' + dQVfk + P^^V^ + i{all - A^W^ + dQV,) = 0. 
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The first coordinate equation gives 

(A.20) 6*fcrO + i{a* - fl + fiq-\u)rl + fiq-Hgi " + alq^^qu^ 

+ blf^ + -fu+ fvQv^QuVj + fUv^dt " qUv^u + a*q-^quy 

= 26jfcr° + i{a* - fi + fiqv^Qu)rl + fiq^^igi - giqv^qu + alq:;^qu)r° 

+ -fu+ fvQv^Qu)rj + fUv^idu - aUv^qu + a*q-\uy = 

For the scalar case, i.e., n = 1, we will denote, for simplicity, a^* = (a^*^, a^*^) by a* = 
(aj, a^). Then, we have 



(A.21) y° 



(A.22) a*=fi- fiq-\. 



(A.23) h*j = -flq^ ^{gi - giq^ ^qu - a*q^ 

= -fiqv^{gi - giqv^qu - {fi - fiqv^qu)qv^qu), 

(A. 24) = (/^g^Hsn - giQv^Qu + alq:;\u) + fvlv^iol - gUv^Qu + a*q~^qu] 



and 

(A.25) y/ 



fiQv ^iai - giqv ^Qu + {ft - fUv ^Qu)qv ^qu) 

+ ft<iv^{9t - aUv^qu + {fi - fiqv^qu)qv^qu] 



1 

-qZ^qu + iqv^{gi - giq^^qu + a*q~^qu)^ 







Li 


H 



Moreover, we let B* be the viscosity matrix [6*^] as in (|1.6p and write 

/ 1 -6** 

(A.26) B* = bl^ 

\b* B* 



□ 
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where b* G M^^^ and B* G Rid-i)x{d~i) _ 



Appendix B. Asymptotic ODE: gap and conjugation lemmas 
Consider a general family of first-order ODE 

(B.l) W'-A(xi,A)W = F 

indexed by a spectral parameter A G C C™, where W G C^, xi G M and ""' denotes 
d/dxi. 

Assumption B.l. 

(ho) Coefficient A(-,A), considered as a function from Q, into C^{xi) is analytic in 
A. Moreover, A(-,A) approaches exponentially to limits A± as xi — > ±00, with uniform 
exponential decay estimates 

(B.2) \{d/dxi)''{A - A±)l < Cie^^l^il/^^^ for xi ^ 0, < A; < K, 

Cj, 6 > 0, on compact subsets ofVL. 

Lemma B.2 (The gap lemma [KS, GZ, ZH]). Consider the homogeneous version F = 0/ 
(jB.ip . under assumption (hO). IfV^{A) is an eigenvector of A- with eigenvalue n{A.), both 
analytic in A, then there exists a solution of (jB.ip of form 

(B.3) W(xi, A) = y (xi, A)e'^(^)^i , 

where V is in xi and locally analytic in A and, for any fixed 9 < 6, satisfies 

(B.4) y(xi,A) = y-(A) + o(e-^"i^ii|y-(A)|), XI <0. 

Lemma B.3. (The conjugation lemma). Given (hO), there exist locally to any given Fq G 
invertible linear transformations P4.(xi,r) = / + 0+(xi, F) and P^{xi,T) = I + 0„(xi, F) 
defined on xi > and xi < 0, respectively, <I>± analytic in F as functions from ^1 to 
C''[0, ±00), such that: 
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(i) For any fixed < 9 < 9 and 0<k<K + l, j>0, 
(B.5) \{d/dAy{d/dxi)''Q±\ < C{j)CiC2e-'^\''^\'^^ for xi ^ 0. 

(a) The change of coordinates W =: P±1^, F =: P±G reduces (jB.ip to 
(B.6) Z' - A±Z = G for xi ^ 0. 

Equivalently, solutions of (|B.ip may be factored as 

(B.7) w=(/ + e±)z±, 

where Z± satisfy the hmiting, constant-coefficient equations (|B.6P and 0± satisfy bounds 

dEl]). 

Example B.4. Consider the linearized equations: 

d 

(B.8) Ut = LU:=- J2i^'U),^ + Qf/ 

(B.9) C/j = L^U := -{A^Uy - ^ i^jA^U + QC/ 

i=2 

Consider a non-homogeneous eigenvalue problem: 
(B.IO) {L^-X)W = f 

Eq. (jB.lOp can 6e expressed in the form: 

d 

W' = -{A^)-\{A^)' + iY,i3A' -Q + ^i)w -{A^r^f 

i=2 

= A(xi,A)VF-F 
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Appendix C. Series expansion of the top eigenvalue of L| 

Consider 

(C.l) L^U = -{A^Uy -iY^^jA^U + QU. 

It can be shown that there exists a unique, analytic eigenvalue 

(C.2) Ao(o = + 7' • I + 1* • {iH) + 

of L| perturbing from the top eigenvalue A = of the operator Lq, with associated analytic 
right and left eigenfunctions 

(C.3) ^{0 = / + ^1 . 1 + iv'l + o(lll') 

and 

(C.4) 7r(0 = 7r° + TT^ • I + e*vr2| + 0(|eT). 

with 99° = tJ' and -k^ = ([n]~\0). 

Lemma C.l. The expansions (fell) . (ICl3]) and (ICH]) /loZd with = a, a G M'^'^ and 
/3 G where —■^'^ = P is positive definite. Here 

(C.5) -d={a2,...,ad) = -{[f^][u]-\...,[f^][u]-') 

Proof. Let ^(1) = 9^0 + Ej^i ^'^j + E,,fc^i ^.^^,4 + ^dlP)- 
The eigenvalue equation L^ip = X^ip leads to 

(C.6) - {A^<^^y - iF\Uy + Q^^ = jjtJ' for j = 2, d. 

Integrating from — cxd to 00 both sides, we have 

_ iF{uy + (0 dg)*«>^d2;i = / jjU'dxi 
-00 J —00 

which yields, by looking at u coordinate, 

(C.8) =7,[C/,] fori = 2,...,d. 
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For u scalar case, we have 

(C.9) a,=i7, = -l^forj = 2,...,d. 

On the other hand, by {'DS) in Assumptions 1.2 and the analytic eigenvalue expansion ()C.2p . 
—7^ = (3 is positive definite. 

□ 
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